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Abstract We report on the optimized production of 38 
a Bose-Einstein condensate of cesium atoms using an 39 
optical trapping approach. Based on an improved trap 40 
loading and evaporation scheme we obtain more than 10 5 41 
atoms in the condensed phase. To test the tunability of 42 
the interaction in the condensate we study the expansion 43 
of the condensate as a function of scattering length. We 44 
further excite strong oscillations of the trapped conden- 45 
sate by rapidly varying the interaction strength. 46 

47 

PACS: 03.75.Kk; 32.80.Pj 

' ■' 49 

1 Introduction 

50 

Much of the present work in the field of quantum gases 51 
relies on optical trapping techniques and on the ability to 52 
tune atomic interactions. Optical approaches have been 53 
recently employed in several atomic Bose-Einstein con- 54 
densation experiments 1,2 3,4 5 and in experiments on 55 
the production of ultracold molecular samples |l)][7||S]0i| 56 
HH] and on molecular Bose-Einstein condensates [1 11112) . 57 
The major advantages in optical traps are the possibility 58 
to trap atoms in any sublevel of the electronic ground 59 
state and the ease to adjust the interaction strength us- 60 
ing magnetically induced Feshbach resonances. 61 
The cesium atom is very attractive for experiments with 
tunable atomic interactions. The lowest internal quan- 62 
turn state of Cs features a unique combination of wide 
and narrow Feshbach resonances which are easily accessi- « 
ble at low magnetic fields ■ This results in a great flex- f>$ 
ibility for tuning the atomic scattering properties. In par- 
ticular, magnetic tuning of the interaction strength has & 
recently allowed the first realization of a Bose-Einstein 
condensate (BEC) with Cs atoms |3] and the realization 69 
of a two-dimensional condensate very close to a dielectric ™ 
surface [5]. The tunability of the atomic interaction can n 
be exploited in experiments where one might wish to ad- n 
just or to dynamically change the mean-field interaction 73 
of the condensate. Also, the Feshbach resonances can be « 
used to produce molecules from an atomic BEC |14ll51 75 

* corresponding author 76 



I51ll0| and to study the transition from an atomic BEC 
to a molecular BEC. In this context, a quantum phase 
transition with an essentially topological character has 
been predicted |15II16| . For such and many other intrigu- 
ing experiments it is desirable to have a large BEC of Cs 
atoms as a starting point. 

In this paper we report on the optimized production of 
an essentially pure Cs BEC in the lowest internal quan- 
tum state with more than 10 5 atoms. Since this state 
cannot be trapped by purely magnetic means, the path 
to condensation relies on a sequence of optical traps. We 
discuss the loading and transfer from one trap to the 
next and give a detailed description of the evaporation 
path and of the resulting condensate. As a demonstra- 
tion for tunability we measure the expansion energy as 
a function of scattering length in time-of-flight experi- 
ments. In particular, we show the ultra-slow expansion of 
the condensate after release from the trap for nearly van- 
ishing scattering length. The release energy corresponds 
to ~ 50 pK. Finally, we present first results when the 
scattering length is suddenly stepped and the conden- 
sate then starts to oscillate freely in the trap. 
2 Cesium scattering properties and Feshbach 
resonances 

Early experiments |17U18j towards condensation of ce- 
sium focused on samples in magnetic traps polarized ei- 
ther in the upper hyperfine ground state F = 4, mag- 
netic sublevel mp = 4, or in the lower hyperfine state 
F = 3, mp = —3. Here, F denotes the total angular mo- 
mentum and mp the magnetic quantum number. The 
spin relaxation rates were measured to be several orders 
of magnitude higher than expected jl9U20U21j . It was 
later understood that this is caused by the dipolar re- 
laxation process induced by the second-order spin-orbit 
interaction |22j . The maximum phase-space density in 
a small sample of Cs atoms was a factor of about four 
away from condensation |23j. 

The problem of the strong inelastic two-body losses can 
be overcome by using the lowest internal state of cesium, 
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Fig. 1 Scattering length as a function of magnetic field for 
the state F = 3, rriF = 3. There is a relatively broad Feshbach 
resonance at 48.0 G due to coupling to a d-wave molecular 
state. The arrows indicate several very narrow resonances at 
11.0, 14.4, 15.0, 19.9 and 53.5 G, which result from coupling 
to g-wave molecular states. The data is taken from |13|. 



F = 3, mp = 3 



. In this state, all inelastic 



two-body processes are endothermic and are thus fully 
suppressed at sufficiently low temperature. This state 
requires optical trapping since it cannot be captured 
in a magnetic trap. Optically trapped atoms can only 
be efficiently evaporated by lowering the total poten- 
tial depth. This process weakens the confinement of the 
trapped sample and thus makes it difficult to achieve 
sufficiently high elastic collision rates for effective evap- 
oration. Hence, adjustability of the collisional properties 
is very helpful for a fast evaporation strategy. 

The success in condensing Cs [2] largely relies on the fact 
that the s-wave scattering length for the F = 3, mp = 3 
state can be tuned to moderate and positive values by 
the application of relatively low dc magnetic fields |13|. 
As Fig. shows, an external magnetic field allows for 
precise tuning of the atomic scattering length a from 
negative to positive values. Positive scattering lengths in 
the range between zero and one thousand ao are attained 
for magnetic fields of a few ten Gauss; ao denotes Bohr's 
radius. In particular, there is a gentle zero-crossing of the 
scattering length near 17 G [21] • Here, the interaction 
of atoms in a BEC is effectively switched off. Several 
narrow higher-order Feshbach resonances ^3] > caused by 
coupling to d- and g-wave molecular states, enable very 
rapid control of the atomic scattering properties. With 
the magnetic field being a free parameter in our optical 
trapping approach, we can take full advantage of this 
tunability of the s-wave scattering length. 

For Cs in the F = 3, mp = 3 ground state the pro- 
cess of three-body recombination is the dominant loss 
and heating mechanism |28j . In a recombination process, 
three atoms collide, two of them form a molecule, and 
the third atom takes away two thirds of the binding en- 



ergy according to energy and momentum conservation. 
The atoms that form the molecule are usually lost, and 
the third atom is either lost or it deposits its share of 
the binding energy in the sample. Heating of the sample 
is the combination of "anti-evaporation" and recombina- 
tion heating |28) . To a good approximation, the three- 
body recombination rate scales with the fourth power 
of the scattering length. Unfortunately, the prefactor in 
this scaling law is measured to be relatively large |28| . 
To minimize this heating, the recombination products 
should be removed quickly from the trap. It is thus im- 
portant to assure that the sample is not operated too 
deep in the hydrodynamic regime and that the evapo- 
ration is efficient in all directions. Arbitrarily increasing 
the scattering length to speed up the forced evaporation 
is therefore not possible without sacrificing cooling effi- 
ciency. Within these limits, tuning the scattering length 
allows for an optimization of the evaporation for given 
trap parameters. For example, for the low initial densi- 
ties in a large reservoir trap the evaporation may be sped 
up by increasing the scattering length. In a later trap- 
ping stage with a higher atomic density the scattering 
length should be reduced to optimize the ratio of good 
to bad collisions. 
3 BEC production 

3. 1 Overview of experimental strategy 

For producing large condensates in optical dipole traps, 
it is necessary to independently optimize both trap load- 
ing and evaporative cooling. For initial loading of as 
many atoms as possible, an optical trap with large vol- 
ume is needed which, in view of limited laser power, im- 
plies a shallow trapping potential. For subsequent forced 
evaporative cooling, however, high densities and fast elas- 
tic collisions require much tighter confinement. These 
two requirements in general demand dynamical changes 
of the trapping potential. A possible way to implement 
this is a spatial compression of the optical trap using 
e.g. a zoom-lens system [21] ■ Our approach is based on 
an alternative way where a sequence of optical trapping 
schemes is used to provide optimized loading together 
with optimized evaporative cooling. 
We first use a shallow, large volume C02-laser trap as a 
"reservoir" for collecting the atoms before forced evap- 
orative cooling is implemented in a tighter trap. The 
reservoir trap can be efficiently loaded with atoms that 
are precooled by Raman-sideband cooling 30 . This ap- 
proach allows collection of atoms at moderate densities 
with little loss from three-body collisions and with neg- 
ligible heating from either photon scattering or trap vi- 
brations. It serves as a good starting point for the final 
transfer into a tighter optical trap. The tighter trap is 
adiabatically increased and adds a "dimple" to the trap- 
ping potential of the reservoir. Collisional loading of this 
dimple already yields a significant enhancement of the 
local number and phase-space density After turn- 
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ing off the reservoir trap excellent conditions for further 57 
forced evaporative cooling are obtained. ss 
The different trap stages of optical trapping used in our 59 
experiments are illustrated in Fig. EI An overview of the eo 
evolution of phase-space density and particle number for 
the various trapping stages is shown in Fig. |3 
The use of relatively weak optical trapping necessitates 62 
the implementation of magnetic "levitation" where a 63 
magnetic field gradient along the vertical direction com- 64 
pensates for the gravitational force. This levitation turns 65 
out to be very useful in two ways: First, in the limit of 66 
very weak optical trapping only one spin state is held 67 
in the trap. This assures perfect spin polarization of the 68 
sample 1 . Further, efficient evaporation can be performed 69 
without the effect of gravitational sag in the trap. The dc 70 
magnetic field offset remains a free parameter for flexible 71 
tuning of the scattering length. 
3.2 Laser cooling 

74 

The initial collection and cooling of Cs atoms is achieved 75 
by conventional techniques. In a stainless steel vacuum 76 
chamber j32j atoms are loaded into a magneto-optical 77 
trap (MOT) from a Zeeman slowed atomic beam with 78 
up to 3 x 10 8 atoms after about 6 s. The MOT is op- 79 
erated on the 6 2 Si/2,-F = 4 to 6 2 P^/2,F' — 5 transi- so 
tion. The ultra-high vacuum of less than 1 x 10 -11 mbar si 
gives 200 s for the 1/e-lifetime of the MOT. The MOT 82 
light is derived from a high power laser diode 2 referenced 83 
via beat-lock to a grating-stabilized master diode laser. 84 
Standard absorption imaging is used to determine par- 85 
ticle numbers and temperatures. se 
We compress the atomic cloud by ramping up the mag- 87 
netic field gradient in the MOT by a factor of 5 to 33 ss 
G/cm within 40 ms. Simultaneously we linearly change 89 
the detuning of the MOT laser from around 10 MHz 90 
to 30 MHz. At the end of the ramp, we switch off the 91 
MOT light and the magnetic field gradient. To cool the 92 
compressed cloud, we then apply degenerate Raman- 93 
sideband cooling [301 i n an optical lattice to further cool 94 
and to polarize the atoms in the desired F — 3, mp = 3 
state. We have adapted the technique as described in 33 
to our setup. This cooling scheme is particularly suited 
for polarizing atoms in the F = 3, mp = 3 state be- 
cause this is a dark state for which photon scattering is 98 
suppressed. Four laser beams derived from an injection 
locked slave laser resonant with the F — 4 to F' = 4 100 
transition produce a three-dimensional optical lattice, 
drive Raman-sideband transitions and repump out of 
the F — 4 ground state manifold. The total power of 103 
all the four beams is 65 mW and their l/e 2 -beam radii 10 " 
are about 1 mm. The oscillation frequency in the lat- 
tice is on the order of 100 kHz. A small magnetic field 106 
offset of several hundred mG is applied to induce the 
Raman-sideband cooling. We succeed in polarizing 90% 108 
of the atoms. The ensemble is then adiabatically released 109 
from the lattice after 6 ms of cooling time. If the atomic 

cloud is released into free space, the temperature of the 111 
~ 112 

1 This Stern-Gerlach separation technique also allows for 
radio-freouenov evaporation alone 1 the vertical direction. Al- 



ensemble with up to 4 x 10 7 atoms is about 0.7 /iK. For 
our typical atomic densities this corresponds to a phase 
space density of 1 x 10 . 
3.3 Reservoir trap 

We generate the large reservoir trap by horizontally cross- 
ing two C02-laser beams Ai and A2 at right angles as 
shown in Fig. \£[&). At the same time we apply a mag- 
netic gradient field in the vertical direction to levitate 
the atoms against gravity. The delivered powers in laser 
beams Ai and A2 are 90 W and 65 W, respectively. 
The light comes from two separate, highly stable lin- 
early polarized single-frequency C02-lasers 3 . Switching 
of the beams is done by external acousto-optical mod- 
ulators 4 (AOMs). Ai is downshifted in frequency by 40 
MHz, whereas A2 is upshiftcd by 40 MHz to prevent any 
interference. To avoid mode-hops the cooling water for 
the lasers needs to be stabilized to better than ±20 mK. 
Still, a slow mode drift changes the power of the lasers 
by a few percent over the time scale of minutes. At the 
crossing point the 1 /e 2 -beam radii of the two lasers are 
(605 ± 35) /im and (690 ± 35) fim. 

The magnetic fields for levitation and for Feshbach tun- 
ing are generated by two pairs of coils aligned with their 
axes parallel to the vertical direction. One pair in anti- 
Helmholtz configuration produces the vertical magnetic 
field gradient near 31.3 G/cm to levitate the atoms in the 
F = 3,mp = 3 state. Another pair in Helmholtz config- 
uration provides a variable bias field Bq of up to 200 G. 
The combined field results in a weak outward directed 
force F(p) — ma 2 p depending on the horizontal distance 
p from the vertical symmetry axis. For perfect levitation 
of our atoms the constant a = g\Jmj (3pbBq) describes 
the curvature of the parabolic anti-trapping potential. 
The levitation field thus slightly reduces the trap depth 
along the horizontal direction. Here, 777 is the mass of 
Cs, g is the gravitational acceleration, and /is is Bohr's 
magneton. At B = 17 G we have a = 2n x 3.4 Hz. The 
horizontal trap frequencies u) x ,y are reduced according 

to uj' x y — \Ju 2 y — a 2 . This is usually a very small effect 

for all but the lowest trap frequencies. Note that levi- 
tation also affects the horizontal motion of free atoms 
after the optical trap is shut off. The horizontal motion 
follows p(t) — po cosh (at) + a~ 1 vo sinh (at) for initial 
position po and initial velocity vq- The vertical motion 
is not affected. 

We excite vertical trap oscillations by briefly changing 
the vertical magnetic field gradient and hence tilting the 
trap. For exciting horizontal trap oscillations we shift the 
equilibrium position of the atoms by adding a horizon- 
tal magnetic field component. In both cases we monitor 
the center-of-mass oscillation of the atomic cloud after 
50 ms time-of-flight. The geometrically averaged trap 
frequency v is calculated to be (12.6 ± 1.5) Hz which 
is in good agreement with the experimental value of 
(13.2 ± 0.2) Hz. Together with the levitation and the 
magnetic bias field the lasers provide an effective trap 

3 Coherent-DEOS GEM-100L 
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Fig. 2 Illustration of the various stages of trap loading and evaporative cooling as seen from above, (a) Plain evaporation in a 
crossed C02-laser trap generated by beams Ai and A2 at a scattering length of a = 1215 ao. (b) 1.5 s of ramping and collisional 
loading into a crossed 1064-nm fibre laser trap generated by beams Bi and B2 with a final scattering length a — 210 ao. (c) 
Forced evaporative cooling after switching off C02-laser beam A2. The power of all remaining lasers is ramped down, and the 
power in C02-laser beam Ai is reduced to zero, (d) Final configuration of the crossed 1064-nm trap. Imaging is done in the 
horizontal plane at an angle of 30° with respect to the long axis of the cigar-shaped atomic cloud. 



depth of about ks x 7/iK. This trap depth is given by 39 
the weaker of the two C02-lasers as the atoms can escape « 
along the direction of the stronger beam. 41 
For transfer of the precooled atoms into the reservoir 42 
trap, we leave the light of the two C02-lasers on dur- 43 
ing the entire pre-cooling phase. This is because the 44 
C02-lasers show strong variations in beam pointing and 45 
beam shape as a function of radio-frequency power to 46 
the AOMs. We have checked that the small light shift 47 
introduced by the lasers does not affect the initial load- 4s 
ing and cooling efficiency. The reservoir trap is then acti- 49 
vated by ramping up the magnetic field and its gradient. 50 
The 1/e-rise time of the magnetic fields is limited to 1.5 51 
ms because of eddy currents in the stainless steel cham- 52 
ber. We therefore do not expect the transfer to be fully 53 
adiabatic. 54 
We find that the atoms are heated to about 2.2 /zK by 55 
the transfer into the reservoir trap. A clear measurement 56 
on the trapped sample is only possible after about 50 ms 57 
since the system initially is not in thermal equilibrium ss 
and since the untrapped atoms need to disappear from 59 
the field of view. We largely attribute the heating to im- 60 
perfect phase space matching. In fact, the atomic cloud ei 
after Raman-sideband cooling to 0.7 /iK has a 1/e-radius 62 
of ~ 350 /im. In comparison, an equilibrium distribution 63 
in the reservoir trap at 0.7 /iK would have a 1/e-radius 64 
of ~ 100 /im. Potential energy is thus gained which is 65 
then turned into kinetic energy, effectively heating the 
cloud of atoms. Subsequently, the hot atoms evaporate 66 
out of the trap. For this phase of plain evaporation we 
set the magnetic bias field to 73.5 G. The scattering len- 68 
gth is then 1215 ao. The temperature is reduced to less 69 
than 1 /xK within 10 s. After this time, we measure more 70 
than 4 x 10 6 atoms, corresponding to a peak phase space 71 
density of 2 x 10 -3 . 
3.4 Dimple trap 

74 

We proceed with loading of the dimple trap after 2 s of 
plain evaporation in the reservoir trap. At this point the 



atom number is 7.8 x 10 6 and the phase space density is 
1.7 x 10~ 3 (see Fig. EJ. The dimple trap is generated by 
horizontally intersecting one tightly focused laser beam 
Bi with 34-jUm waist and another less focused beam B2 
with 260-/im waist at right angles, rotated by 30° in the 
horizontal plane with respect to the COvlaser beams 
as shown in Fig. Eld). This is different from our ear- 
lier work 0] where we have used C02-laser beam A2 for 
axial confinement. We introduce the B2 beam because 
some weak back reflections of the COvlaser beams led 
to a slight undesirable corrugation of the optical poten- 
tial. This complicated the quantitative analysis of the 
BEC. Beams Bi and B2 are derived from a broadband 
fiber laser 5 at 1064 nm. The powers in these beams are 
ramped up within 1.5 s to a maximum power of 70 mW 
for Bi and 270 mW for B2. The trapping in the dimple is 
now briefly done by all four laser beams with Bi provid- 
ing most of the radial and Ai most of the axial confine- 
ment. After switching off beam A2 we measure the radial 
and axial trap frequencies in the dimple to (221.2 ± 1.6) 
Hz and (14.2 ± 0.1) Hz, respectively. During the ramp- 
ing up phase of Bi and B2 we reduce the magnetic field 
offset to 23 G and thus the scattering length to 300 ao 
in order to reduce losses from three-body recombination 
|28| . The trap now contains about 1.7 x 10 6 atoms at a 
peak phase space density of approximately 0.13. 
3.5 Forced evaporation towards BEC 

We start forced evaporative cooling by ramping down 
the power in all three remaining beams. Simultaneously 
we remove the reservoir by switching off the C02-laser 
A2 that is not responsible for axial confinement. To as- 
sure a well-defined ramp over a large intensity range we 
control the light power of the near-infrared beam Bi by 
means of a logarithmic photodiode and a servo loop. The 
power in C02-laser beam Ai is ramped to zero within 
5.5 s so that B2 at the end of evaporation exclusively as- 
sures axial confinement. The change in beam pointing for 
A2 does not affect the evaporation. For Bi we approxi- 

5 IPG Laser PYL-10 
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atom number 

Fig. 3 Peak phase space density as function of atom num- 
ber. The path of evaporation proceeds from right to left. The 
triangle shows the atomic ensemble immediately after lattice 
cooling. The open circles show the ensemble in the reservoir 
trap after 0.08, 0.22, 0.64, and 2.0 s. The filled circles corre- 
spond to the sample in dimple trap right after loading and 
after 1.5 s of evaporation. The phase transition occurs after 
2 s of forced evaporation with ~ 5 x 10 5 atoms left in the 
dimple trap. 



mately follow an exponential ramp over 5.5 s. The power 
in beam B 2 is only slightly reduced. The final power in 
Bi and B2 is 0.5 mW and 220 mW. We find and optimize 
this ramp by extending the ramp in discrete time steps 
of a few hundred milliseconds at the beginning and up 
to one second towards the end of the ramp. 
At each step we search for a maximum in evaporation 
efficiency 7 = log(D' /D) / \og(N/N') as a function of the 
trap depth and scattering length [331. Here, D and D' 
are the phase-space densities at the beginning and end 
of each step, N and AT' denote the respective particle 
numbers. Maximizing 7 at each step results in an overall 
optimization of the evaporation path. We find that a 
magnetic field value of 21 G with scattering length a = 
210 ao is optimal during the forced evaporation phase. 
As can be seen from Fig. [21 the efficiency 7 lies around 
3 during the forced evaporation ramp. We attribute this 
high efficiency to the fact that atoms can escape the trap 
into almost all directions because of the levitation field. 
We observe the phase transition after 2 s of forced evapo- 
rative cooling with about 5 x 10 5 atoms at a temperature 
of (200 ± 10) nK. At this point the power in beams Bi 
and B 2 is 8.7 mW and 250 mW. The duration of the 
ramp is relatively short. Our evaporation proceeds close 
to the hydrodynamic regime. Thus, significant improve- 
ment of the evaporation is not to be expected. 
Further evaporation leaves a cigar-shaped condensate 
with the long axis in the horizontal plane. In Fig. 0] we 
show vertical density profiles of expanding condensates. 
The tunability of the scattering length allows us to ex- 
plore different regimes of expansion. For Fig. ^a) we 
expand the condensate at the creation scattering length 
of 210 ao- This is the usual type of self-similar expansion 
in which the condensate in the Thomas-Fermi regime re- 



o 

-D 
"i_ 

GO 

T3 
>* 
CO 

CD 
"O 

"O 

CD 
N 

E 
o 



1.0 



0.5 



0.0 



o 

' 

Z3 
-Q 
"i_ 

CO 
T3 

>s 
■ H — ■ 

CO 

CD 
"O 

U 
CD 
N 

"CO 

E 
o 



1.0 



0.5 



0.0 





-300 -150 150 

vertical position (\im) 



300 



Fig. 4 Vertical density profiles of Cs condensates after 100 
ms of free expansion in the levitation field. The solid curves 
are fits to the data for the Thomas-Fermi profiles which in- 
clude possible thermal components. For better distinction the 
baseline is dashed, (a) Expansion with no change in scatter- 
ing length. The total number of atoms in the condensate is 
N = 1.1 X 10 5 . (b) Expansion near zero scattering length un- 
der the same conditions reveals a small thermal component 
with a temperature of about 10 nK. 



tains its parabolic shape [35] • For Fig. 0Jb) we step the 
scattering length to zero at the moment of release from 
the trap. The mean-field interaction thus vanishes and 
the rate of expansion is greatly reduced. This exposes a 
small thermal component, for which a bimodal fit reveals 
a temperature of around 10 nK. The critical tempera- 
ture at these trapping conditions is 24 nK, therefore the 
expected condensate fraction agrees well with the mea- 
sured value of 91%. From the fit to the data in Fig. 01 
we obtain that there are up to 1.1 x 10 5 atoms in the 
condensate with a 20% calibration error. The error does 
not come from the fit but from the overall uncertainty in 
determining the atom number. Usually, the error from 
absorption imaging alone is around 50%, but we can cali- 
brate the atom number from measurements on the chem- 
ical potential, see Sec. 14.11 For this particular experiment 
we measure the final trap frequencies to (4.3 ± 0.2) Hz 
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and (21.1 ± 0.2) Hz along the axial and radial direction, 
respectively. We thus infer for the initial Thomas-Fermi 
sizes R TF = (8.7 ± 0.3) fim and R TF = (42.5 ± 1.2) ^m 
along the radial and axial directions at a scattering len- 
gth of a = 210 do- The peak density of the condensate is 
n = (2.1 ±0.1) x 10 13 cm~ 3 . 
4 Tunable quantum gas 

We now test the tunability of the condensate interac- 
tion. We first study the condensate expansion as a func- 
tion of scattering length in two different ways. We 
then specialize to the case when the interaction energy is 
switched off and present improved results on the ultra- 
slow expansion of the condensate in comparison with 
earlier measurements in ^Q. Finally, we excite compres- 
sion oscillations of the trapped condensate by suddenly 
stepping the scattering length to a lower value. 
4-1 Expansion energy as a function of scattering length 

We measure the release energy of the condensate for 
slow and fast changes of the scattering length. When we 
slowly vary the scattering length the wave function of 
the trapped condensate can follow adiabatically and the 
condensate remains in equilibrium. The release energy is 
proportional to the chemical potential of the condensate 
at the given value of the scattering length. The situation 
is different when we rapidly switch the scattering length . 
at the moment of condensate release. The condensate 
then expands from a non-equilibrium state because the 
wave function has not had time to adjust to the change 
in interaction energy. This leads to strong changes for 
the rate of condensate expansion in comparison to the 
expansion from equilibrium. 

We first consider a condensate in the Thomas-Fermi 
regime for which we adiabatically ramp the scattering 
length to a new value. For such a condensate, the re- 
lease energy E re i directly corresponds to the chemical 
potential htf through ^E re i = ^tf G2j, which is given 

by 



hv flhN 



2/5 



l*TF = ^ — « 2/5 - (1) 

2 V a ho J 

Here, v is the geometric average of the trap frequen- 
cies, N is the particle number in the condensate, and • 

= \/h/(m2TTp) is the oscillator length. For the ex- 
periment we produce a condensate with N = 8.5 x 10 4 • 
atoms at a creation scattering length of a c = 210 ao- We 
then slowly ramp the magnetic field to values between 
20 and 35 G, setting the scattering length to a value be- 
tween about 200 and 700 ao. The slow ramping excludes 
values below the Feshbach resonance at 19.9 and above 
the one at 48.0 G because of strong loss 6 . The conden- 
sate is then released from the trap and we measure the 
release energy. The results are shown in Fig. [S] Here we 
assume that the magnetic field strength translates into 
scattering length according to Fig. ^ The data is well 
fit by a function of the form C a 2//5 according to Eq. (JTJ. 
From the fit parameter C we can deduce an independent 

6 A combination of slow ramping and quick jumping at the 
Feshbach resonances would allow access to the full range of 
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Fig. 5 Release energy of the condensate as a function of scat- 
tering length a. The filled circles represent experimental data 
for the case of adiabatic ramping of a trapped condensate. 
The data, corresponding to 2/7 of the chemical potential at a 
given value of the scattering length, is fit by Ca 2//5 . The open 
circles represent data for rapid switching at the moment of 
condensate release. As discussed in the text, the straight line 
is not a fit. It connects the origin with the fitted value of the 
release energy at the creation scattering length. 



estimate of the particle number N = (8.2 ± 1.3) x 10 4 . 
The error is dominated by the error in determining the 
trap frequencies. 

For a sudden change of the scattering length the conden- 
sate wave function has no time to react. For example, 
for an increase of the scattering length the density dis- 
tribution is too narrow in comparison to the equilibrium 
density distribution at the new value of the scattering 
length. The condensate thus expands more rapidly than 
a condensate in equilibrium at this new value. Since the 
mean-field interaction energy of the condensate scales 
linearly with the scattering length for a given density 
profile |35| . we expect a linear behavior of the release 
energy as a function of the final scattering length a. In 
Fig.|S]we thus compare the data for the measured release 
energy to a straight line C a 2 J 5 a/a c given by the origin 
and the fitted value of the release energy at the creation 
scattering length a c — 210 a . We find good agreement 
with the linear dependence. 
4-2 Ultra-slow condensate expansion 

We now study the expansion of the condensate near the 
zero-crossing of the scattering length. At the moment of 
condensate release, we rapidly switch the magnetic field 
from the creation field near 20 G to (17.17±0.05) G, cor- 
responding to a — (3.4±3.0) ao- The error in determining 
the precise magnetic field at the position of condensate 
requires that we choose a slightly positive value of the 
scattering length to assure that no weakly attractive in- 
teractions modify the condensate expansion. The levi- 
tation field remains on, allowing for an extended obser- 
vation period because the atoms then do not fall under 
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80 100 120 140 

time of flight (ms) 



Fig. 6 Expansion of the non-interacting condensate. The 
data points show the horizontal (above) and vertical (be- 22 
low) rms radius of the BEC as a function of expansion time 23 
near the zero crossing of scattering length. Note the different 24 
scales. The fit to the residual vertical expansion reveals a re- 
lease energy of ks x (51 ±3) pK. For the horizontal expansion 
the data is fit by A cosh (at) with a = 2tt x (3.20 ±0.23) Hz. 

25 
26 
27 
28 

gravity. Fig. [5] shows the vertical and horizontal extent 29 
of a BEC with 1.2 x 10 5 atoms as a function of time after 30 
release from the trap. We only show the data after 75 31 
ms of expansion when the optical density of the atomic 32 
cloud is sufficiently reduced to allow for reliable absorp- 33 
tion imaging. The horizontal expansion is dominated by 34 
the magnetic anti-trapping potential which derives from 35 
the presence of the levitation field and which magnifies 36 
the atomic cloud according to the cosine hyperbolicus 37 
function, see Sec. 13.31 The measured rate of expansion 38 
2tt x (3.20 ± 0.23) Hz agrees reasonably well with the 39 
expected rate constant a = 2ir x 3.4 Hz. The vertical ex- 40 
pansion corresponds to a release energy of fee x (51 ± 3) 41 
pK. Note that this is much lower than the kinetic energy 42 
of the ground state hw r /A — ks x 253 pK given by a ra- 43 
dial trap frequency of u> r = 2tt x 21.1 Hz. It is remarkable 44 
that the release energy is less than the zero-point energy 45 
of the ground state. Since the spatial extent of the con- 46 
densate is much larger than the size of the ground state 47 
wave function of the harmonic oscillator, the momentum 48 
spread, limited by the uncertainty of the wave function 49 
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30 t , , ^ 

50 100 150 

hold time (ms) 

Fig. 7 Condensate oscillations after rapid switching of the 
scattering length. The filled circles show the vertical rms ra- 
dius of an expanding BEC with 7 x 10 4 atoms after 80 ms 
of free expansion as a function of hold time in the trap. The 
scattering length has been switched rapidly from 363 ao to 
25 ao. The solid curve is a fit to the data giving an oscillation 
frequency of (58.3 ± 0.2) Hz. We independently measure the 
radial trap frequency to (28 ± 1) Hz. 

of the initial condensate, is lower than that of the ground 
state. 

4-3 Condensate oscillations 

By rapidly ramping the scattering length it is possible 
to excite oscillations of the condensate in the trap |37| . 
In fact, in the limit of a cigar shaped condensate one ex- 
pects radial "compression" or "expansion oscillations" at 
twice the trap frequency. Compression oscillations can 
be seen in Fig. where we plot the vertical radius of 
the released condensate as a function of hold time in 
the trap. To excite the oscillation we step the scatter- 
ing length from a value of a — 363 ao (B — 24.4 G) to 
a = 25 ao (B — 17.6 G) at time to- The condensate is 
then allowed to oscillate in the trap for a variable hold 
time th at the final value of the scattering length. We 
release the condensate at time to + th and take an im- 
age after 80 ms of free expansion. We fit the data by a 
sinusoidal function. The measured compression oscilla- 
tion frequency of (58.3 ± 0.2) Hz agrees well with twice 
the radial trap frequency of 2 x (28 ± 1) Hz at the given 
trapping power. To account for the damping we have to 
introduce an exponential decay of the amplitude and of 
the offset value. The damping of the amplitude has a 
time constant of 126 ms. We have not yet identified the 
origin of this damping. Possibly the BEC samples differ- 
ent trapping frequencies due to the large amplitude of 
the oscillation, which would lead to an apparent damp- 
ing. Also, damping might be caused by the interaction 
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with a residual thermal cloud or by parametric processes 58 

[SHI - 

5 Conclusion 60 

61 

We have shown that essentially pure Cs condensates can 62 
be produced with more than 10 5 atoms. In our optical 
trap it is possible to flexibly change the atomic scatter- 65 
ing properties. The atomic condensate can now be used 
as the starting point for experiments where a tuning and 67 
ramping of the scattering properties can be exploited. It 6S 
will be interesting to study the case of a non-interacting 69 
condensate at the zero-crossing of the scattering length. ™ 
Such a condensate might be used in atom interferometers 71 
where one wishes to suppress any mean- field effects jSH] • 72 
On the other hand, tuning to large values of the scatter- 73 
ing length might allow the investigation of effects beyond 74 
the mean- field approximation |35) . Also, modulation of 
the scattering length could be used as an alternative tool ?? 
to probe the excitation spectrum of the condensate. Fi- 78 
nally, ultracold Cs2 molecules can be created by ramping 79 
across one of the Feshbach resonances [B] and the tran- 8 o 
sition from an atomic to a molecular condensate could si 
then be studied. 82 
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